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I  .      INTTRODUCT  ION 
A.      OVERVIEW 

In  the  design  process  of  composite  material  structures  it 
is  often  noticed  that  the  best  characteristics  of  the  material 
are  not  used.  The  normal  procedure  is  to  tailor  the  material 
to  approach  the  properties  of  isotropic  materials,  avoiding 
the  "undesirable"  couplings  between  extension  and  twist  for 
example . 

That  traditional  approach  worked  for  helicopter  design. 
However,  in  the  tilt  rotor,  the  idea  of  changing  the  twist 
angle  of  t  lie  blade  to  optimize  for  either  hover  or  forward 
flight  came  into  the  scenario. 

This  study  will  use  a  blade  spar  to  study  the  sensitivity 
of  the  extension-twist  coupling  due  to  changes  in  RPM ,  layer 
angles  definition,  changes  in  laminate  orientation  and  also 
the  effects  of  each  of  these  configurations  in  the  first  four 
natural  vibration  modes  of  the  blade. 

In  the  analysis,  a  program  named  CASA/GIFTS  (Computer 
Aided  Structural  Analysis/Graphical  Interactive  Finite  Element 
Total  System)  is  adopted.  The  program  makes  it  easier-  to 
visualize  the  effects  of  the  coupling  in  the  structuri  and  is 
--  en  to  !  an  efi  icient  design  tool.  The  results  obtained 
through    h    I  Li   t<  element   procedure   for   ea  :    Laminate 
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configuration,  in  a  total  of  six  models  will  be  presented  in 
3-D  graphs  that  show  the  trends  for  each  model  for  the  layer 
angle  (top  layer  as  a  reference),  the  RPM  and  the  twist  angle. 
These  graphs  can  then  be  used  as  a  preliminary  orientation 
by  the  designer  to  achieve  the  optimal  or  near  optimal 
combination  of  laminate  configuration  layer  angle  within  the 
laminate  that  fit  the  specific  requirements  for  the  particular 
blade  in  focus . 

B.      LITERATURE   REVIEW 


This  section  reviews  briefly  some  pertinent  publications 
listed  as  References  in  the  thesis  report.  Other  related 
literature  found  during  the  scope  of  this  research,  is  listed 
under  Bibliography.  The  topics  reviewed  are  classified  under 
the  following  categories: 

Finite  Element  Theory. 

Composite  Material  Theory. 

Vibrations . 

Helicopter  Theory. 

Design  Applications  (Modeling). 

Auxiliary  Software. 

In   the   Finite   Element   Theory,   relevant   to   present 

research,   some   basic   references   are:   Batozf 19801   that 

describes   the   finite   element   procedures   in   engineering 

analysis  using  Discrete  Kirchoff  Theory  ( DKT  )  ,  Cookf 19811  on 


basic  concepts  and  applications  of  finite  element  analysis, 
Craig[1981]  describing  computer  methods  in  structural  dynamics 
and  fCASA/GIFTS,  1987]  that  is  the  reference  manual  for  the 
finite  element  code  used  for  the  present  ( static  and  dynamic  ) 
work.  Its  graphical  capabilities  are  very  convenient  to  the 
analyst  in  a  manner  that  simplifies  the  visualization  of  the 
effects  of  loads  and  couplings  in  the  structure,  making 
difficult  concepts  more  accessible  to  understanding. 

Tsai  and  Pagano[1968]  establish  a  notation  in  which  the 
composite  lamina  properties  are  invariant  with  respect  to  the 
axis  of  rotation.  This  approach  provides  a  very  useful  way  to 
compare  various  material  systems;  the  laminate  theory  is  well 
documented,  see  for  example  Vinsonfl987]  and  JonesT 19751  both 
offer  a  good  source  on  the  mechanics  and  behavior  of  composite 
materials . 

Yntema[ 19551  developed  a  very  useful  tool  to  estimate 
bending  frequencies  of  rotating  beams  in  which  the  stiffness 
effect  of  RPM  is  shown  for  several  blades  (beam) 
configurations.  Schilhansl \ 1958  ]  uses  another  approach  to 
estimate  bending  frequencies  of  rotating  cantilever  beams, 
introducing  the  effect  of  the  angle  made  by  the  minor  axis  of 
inertia  with  the  direction  of  the  tangential  circular 
velocity.  Woodf 19651  shows  a  parametric  investigation  of 
aerodynamic  and  aeroelastic  characteristics  of  articulated  and 
hingeless  rotor  systems.  Pr i tchardf 1 988 1  presents  an  optimal 
placement  of  tuning  masses  to  reduce  vibration  levels  in 


helicopter  blades.  Kottapall i [ 1 983 1  approaches  the  vibration 
reduction  problem  by  modifying  the  blade  torsional  response. 
This  is  an  innovative  approach  in  that  he  modifies  the 
aerodynamic  response  instead  of  adding  mass,  which  is  the 
usual  conventional  way  to  solve  the  vibration  reduction 
problem. 

The  helicopter  theory  in  general  is  found  comprehensively 
in  Johnsonf 1 980  ]  and  is  complemented  in  this  research  by 
Prouty[1986]  in  the  understanding  design  trends. 

Hoskin[1986]  is  a  general  informative  book  in  the 
aeronautical  utilization  of  composite  materials,  and  in 
Lake[1988]  there  is  a  preliminary  investigation  of  finite 
element  modeling  of  composite  material  rotor  blades. 
Nixonfl987]  addresses  the  extension-twist  coupling  of 
composite  circular  tubes  applied  to  tilt  rotor  design.  In 
McVeigh  [1983]  the  aerodynamic  design  of  XV-15  advanced 
composite  tilt  rotor  blade  requirements  are  reviewed.  Finally 
Hodges[1987]  presents  a  comparison  of  composite  material  rotor 
blade  models,  using  two  different  methods  of  analysis:  this 
model  was  chosen  as  a  prime  reference  to  build  the  basic  model 
used  throughout  this  work. 

Several  software  tools  made  it  possible  to  manipulate  the 
immense  volume  of  data  transforming,  it  into  a  suitable  and 
compact  presentation  form,  enabling  the  conclusions  to  be 
obtained.  Kelly [1988  1  was  the  basic  reference  to  use  Word 
Perfect  5.0  in  the  word  processing  and  graphical  editing  jobs. 


Holtf 19881  is  a  reference  that  helps  to  manage  the  data  into 
Lotus  1-2-3  files  and  to  get  the  input  graph  files  to  Word 
Perfect  and  to  SURFERf 1987 1 .  The  SURFER  program  was  used  to 
obtain  the  3-D  carpet  plots  of  extension-twist  coupling 
analysis . 

C.      THESIS   OUTLINE 

This  thesis  report  is  divided  into  three  areas,  often 
found  during  a  design  process. 

It  begins  by  modeling  a  rotor  blade  with  laminate 
composite  material.  A  numerical  solution  is  obtained, 
exploiting  the  extension-twist  coupling  inherent  when 
asymmetries  are  present  in  the  laminate  construction. 
Discrepancies  encountered  in  the  results  compared  with  other 
rotor  models  are  reported.  These  differences  raised  questions 
fostering  the  research  towards  the  possibilities  of  different 
rotor  constructions,  starting  with  the  same  basic  laminate. 

The  different  rotor  configurations  obtained  are  then 
subject  to  both  static  and  dynamic  analysis.  The  static 
analysis  reveals  the  different  responses  of  these  models, 
while  the  dynamic  behavior  remains  almost  invariant,  yielding 
only  small  variations  in  the  first  four  fundamental  vibration 
modes . 

Analysis  of  the  results  is  done,  and  it  is  believed  that 
a   new   insight   into   finite   element   analysis   of   composite 


laminates  is  provided.  This  new  insight  is  better  noticed  when 
there  are  asymmetries  present  within  the  laminate 
construction . 


1 1  . 


THEORETICAL   FORMU L AT I ON 


A  . 


BASIC   CONS I DERATIONS 


To  understand  how  composite  material  works  and   its 
particular   behavior   that   is   useful   in   "designing"   the 
material,  a  macroscopic  approach  is  resorted  to. 

Composite  means  that  two  or  more  materials  are  combined 
and  usually  preserve  their  best  qualities,  and,  often  exhibit 
properties  that  neither  one  possesses  by  itself. 

Some  of  the  properties  that  can  be  improved  by  making  a 
composite  material  are: 


strength 
stiffness 

corrosion  resistance 
wear  resistance 
life  in  fatigue 
thermal  insulation 
thermal  conductivity 
acoustical  insulation 


It  may  be  noted  that  not  all  these  properties  are  improved 
at  the  same  time,  and  usually  it  may  not  be  required. 


An  Orthotropic1  material  has  properties  which  are 
different  in  three  mutually  perpendicular  directions  at  a 
point,  in  other  words  has  three  different  planes  of  symmetry. 

A  lamina  is  the  basic  form  of  any  laminated  fiber- 
reinforced  composite,  it  consists  of  a  plane  mounting  of 
fibers  imbedded  in  a  matrix. 

In  the  macro-mechanical  behavior  of  the  lamina,  only 
averaged  apparent  mechanical  properties,  or  better,  the 
stress-strain  relations  for  an  orthotropic  material  under 
plane  stress  conditions,  transformed  to  directions  not  aligned 
with  the  principal  directions  of  the  lamina,  are  considered. 

This  describes  the  behavior  of  the  orthotropic  material 
whose  laminae  have  arbitrary  directions,  other  than  the 
natural  geometric  directions  of  the  structure  to  be  analyzed. 

B.      STRESS— STRAIN   RELATIONS 

The  stress-strain  relations  or  generalized  Hooke ' s  law  for 
anisotropic  material  will  be  reviewed  briefly,  and  the  basic 
restriction  of  linear  elastic  behavior  assumption  applies. 

The  following  derivations  can  be  found  in  more  detail  in 
Vinson! 1987]  . 

The  elasticity  tensor  has  nine  independent  elasticity 
constants,  because  of  symmetry  (0^=0^)  and  the  three  mutually 


1-    Orthogonally  anisotropic. 
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orthogonal    planes    of    symmetry    ( C16=C26=C36  =  C45  =  0 )  ; 
that    is , 


(1) 


Cii 

c12 

Cl3 

0 

0 

0 

c21 

c22 

C23 

0 

0 

0 

i  J 

C31 

C32 

^33 

0 

0 

0 

0 

0 

0 

c44 

0 

0 

0 

0 

0 

0 

C55 

0 

0 

0 

0 

0 

0 

c 

66 


and 


(2) 


o,    =    C±*   x    e<       (1  =  1,2,3,4,5,6;    j  =  l ,  2  ,  3  ,  4  ,  5  ,  6  ); 


ij 


In  these  relations  the  hygrothermal  and  thermal  effects 
are   not  included  for  the  sake  of  simplicity  of  this  report. 

Through  a  force  equilibrium  study,  analogous  to  the  Mohr's 
circle  analysis  [Jones,  1975]  the  stresses  and  strains  in  the 


2The  material  directions  1  to  6  are  equivalent  to 
cartesian  directions:  x  ,  y , z , yz , xz ,xy  respectively  in  material 
ax  i  s  . 


material  directions  can  be  related  to  arbitrary3  directions 
using  the  following  relations: 


(3) 


" 

°x 

°y 

az 

ayz 

°xz 

°xy 

=  [T]"1 


and 


(4) 


yz 


■xy 


=  m 


-i 


e2 

e3 


The  transformation  matrix  is  constructed  by  defining 
m=cos  0,  n  =  sin  0;  where  0  is  the  ply  orientation  angle, 
defined  with  respect  to  material  and  structural  axes,  assuming 
positive  sign  in  counterclockwise  direction  (Figure  1). 


3  More  convenient,  xyz  axis  related  to  the  structure 
geometry  for  example. 
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Figure  1  Positive  rotation  of  principal  material  axes  from 
arbitrary  xy  axes    [Jones , 1975] . 


The  inverse  of  transformation  matrix  [T]  is  expressed  in 
terms  of  direction  cosines  m  and  n  as: 


(5) 


[T] 


-1 


m2 
n 

n2 
m 

0 

0 

0 

-2mn 

0 

0 

0 

+  2mn 

0 

0 

1 

0 

0 

0 

0 

0 

0 

m 

-n 

0 

0 

0 

0 

n 

m 

(m    -    n    ) 

mn 

0 

0 

0 

0 
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For  a  lamina  of  composite  material  the  modern  notation 
uses  Qij  for  the  quantities  in  the  material  stiffness  matrix 
and  are  defined  as: 


Q11=E11(1-V23V32)/5'  Q22=E22(1-V31V13)/6 

Q33=E33(1-V12V21)/6'  Q44=  G23'    Q55=  G13 

(6)  Q12=E11(V21+V31V23)/6  =  E22(v12+V32V13)/6 

Q13=E11(V31+V21V32)/6  =  E22(V13+V12V23)/6 

Q23  =  E22(V32  +  V12V31)/6  =  E33 ( V23  +  V21V13  >  /& 


fi=    —  v   v   —  v   v   —  v   v   —  2v   v   v 

12  21    23  32    31  13     21  32  13 


If  the   lamina  has   same   properties   in  both  2   and   3 

directions,    then  vi2  =  vi3>   v3i  =  v32  =  0>   ^12=(^13»   ^22  =  ^33  an<^  a 
simplification  in  equations  (6)  lead  to  equations  (7). 


Qn=En/(l  -  v12v21),      Q22  =  E22/(1  "  V12V21 


(7)  Q12=    v21En/(l     -    v12v21)=    v12E22    /(I     -    v12v21) 


Q66=    G12 
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To  compare  various  material  systems  in  the  design  of  a 
composite  structure,  Tsai  and  Pagano[1969]  arrived  at  some 
invariants  with  respect  to  the  axis  of  rotation: 

V  I  (3Q11+3Q22+2Q12+4Q66) 

V  £  <Q11-  Q22> 

(8)  U3=  |  (Qn+  Q22-  2Q12-  4Q66) 

V  I  ,Q11+  Q22+  6Q12-  4Q66> 


V  I  (Q11  +  Q22-  2Q12+  4Q66» 


These  invariants  are  representative  of  the  lamina 
properties.  The  [Q]  matrix  may  be  transformed  to  x-y  axes  from 
1-2  axes  by  the  usual  coordinate  transformation. 

(9)  [Q]  =  [T]"1  [Q][T] 

In  terms  of  the  invariants  Us,  the  fQ]  matrix  elements 
are  written  as : 


Qll=    Ul+    u2cos(2e)    +    U3cos(49) 


Q22=    l^-    U2cos(2G)    +    U3cos(46) 


Q12=    U4-    U3cos(46) 


(10) 


Qcc=    un"    U0cos(46) 
66  5  3 

Q16=    |    U2sin(2G)    +    U3sin(46) 
Q26=    |    LT2sin(26)    -    U3sin(40) 
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After  the  definitions  in  (10),  the  stress-strain  relations 
for  the  kth  lamina  of  a  N  ply  laminate  (Figure  2)  are  given  by: 


(11) 


ryz 

xz 

r 
xy 


=  [Q] 


Jk 


x 


.yz 

ixz 

:xy 


-»k 


Where  Qi;j  are  defined  earlier. 


Figure  2  Laminate  construction,  [ 90°/+45o/0o/-45o/90° ] 
characterized  in  this  case  to  show  the  nomenclature  for 
the    lay-up. 


14 


C.      LAMINATE   ANALYSIS 

When  two  or  more  laminae  are  bonded  together  to  obtain  an 
integral  structural  component,  the  laminae  principal  material 
directions  are  established  to  produce  a  component  that  can 
resist  loads  in  several  directions  (Figure  2). 

The  following  procedures  enable  to  obtain  the  stiffness 
matrix  of  such  a  composite  material  from  the  lamina 
properties,  using  arbitrary  combinations  of  materials  and 
layer  angles. 

For  the  kth  lamina  of  a  N  layer  laminate,  equation  (11) 
represents  the  stress-strain  relations. 

In  order  to  obtain  the  strain-displacement  relations,  the 
displacements  are  defined  as  follows: 


u(x,y,z)    =    u0(x,y)    +    za(x,y) 
(12)         v(x,y,z)    =    v0(x,y)    +    zR(x,y) 
w(x,y,z)    =    w(x,y) 


where  u,  v   and  w  are  displacements  in  x,  y  and  z  directions 
respect  ively 
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and' 


(  13  ) 


6w  6w 

a  =  -  —    and    ft  =  -  — 

8x  6y 


are  the  negative  of  the  slope,  related  to  the  rotations 


The  mid-surface  strains,  then  are  given  by 


x 


6u 
o    6x 


o   ' 


5vQ  i 

6y 


(  14 


1  f  6u       6v 

u  oy       ox   J 


while  the  curvatures  are  given  by 


K 


6a  56 

6x     '  K  y=    6y 


(15  ) 


xy 


1  I"   6a      Ml 

2  6y      6x 


4   In   the 
derivative . 


following   expressions   6   means   partial 
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Defining  h  as  laminate  thickness,  hk  (Figure  3)  as  the 
vectorial  distance  from  the  panel  mid-plane,  N,  M,  Q  as  stress 
resultants  (Figure  4),  stress  couples  and  shear  resultants 
respectively.  Relations  can  be  established  relating  these 
resultants  to  the  strains  and  curvatures.  This  matrix  is 
henceforth  referred  as  the  ABD  matrix. 

This  nomenclature  is  adopted  for  clarity  and  simplicity 
as  may  be  noted  from  the  expressions  that  follow. 

The  stress  resultants  are  expressed  by  integrating 
stresses  across  the  thickness  of  the  plate  and  is  given  by: 
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xy 
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xy 
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and  the  moment  resultants  take  the  following  form 


(17) 
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Figure  3    Nomenclature        for        the        Stacking        Sequence 
[Vinson,  1987]. 


Figure  4  Positive    directions     for    Stress    Resultants    and 
Stress    Couples    for   a    Plate    [Vinson , 1987 ] . 
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Relating  resultants  with  the  inplane  strains  and 
curvatures  through  expressions  (11)  through  (14),  matrix  form 
is  given  by: 


(18) 


[N]  =  [A].[eJ  +  [B].[K] 


[Ml  =  [B].[€0]  +  [D].[K1 


In  expression  (18),  the  matrices  A,  B  and  D  are  called 
extensional  stiffness  matrix,  bending-stretching  coupling 
matrix  and  flexural  stiffness  matrix  respectively,  where: 


(19) 


Aij  =  2k=l  (  *ij>k  [  hk  "  hk-l] 

B±j  =  i  *  -Li(  ho  \  r  hk  -  hk-i] 

DiJ    =   1   *  2k=l(    *ij    >k   [    hk  "   hk-l] 


[i,j=l,2,6] 
[i,j=l,2,6] 
[i,j=l,2,6] 


Assuming  that  transverse  shear  stresses  are  distributed 
parabol ically  (even  though  there  exists  discontinuities  at 
the  laminae  interfaces ) ,  the  shear  resultants  may  be  written 
as : 


(  20 


0   =  2(  A,K6    +  A  .-6 
x        5  5  x  z     4  5  vz 


Q   =  2(  A,_€    +  A,  .6   ) 
y        4  o  x z     44  yz 
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where 


(21  ) 


\r  isk=i  '  «ij»k 


V  hk-r  I  *  ,hk"  hk-i'*  h 

n 


i,j  =  4,5  only 
The  stresses  and  moments  resultants  can  be  written  in  a 
compact  form  that  is  very  convenient  to  see  the  overall 
relations  within  a  composite  material  laminate,  specifically 
the  ABD  matrix.  The  ABD  matrix  is  obtained  as  a  combination 
of  the  expressions  in  (18).  Expressions  (18)  and  (22)  are 
equivalent  and  in  expression  (22)  the  ABD  matrix  is  defined 
explicitly. 
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EFFECTS   OF   FLY   ANGLE   O 


To  get  an  insight  of  the  overall  influence  of  the  ply 
angle  6  on  the  laminate  stiffness,  it  is  useful  to  observe 
the  effects  of  G  for  a  given  lamina.  As  each  element  of  the 
ABD  matrix,  Ai1?  Bijt  and  T>i),  is  a  function  of  the  material 
stiffness  coefficient  QA1,  the  layer  thickness  tk  and  mid-plane 


2  0 


vectorial  distance  hk,  the  Qij's  represent  the  stiffness 
variations  due  to  ply  angle  changes  in  a  given  lamina. 

The  property  of  the  composite  material  that  reveals  its 
advantages  over  metallic  or  plastic  is  that  it  can  be  tailored 
to  fit  a  given  structural  shape  and  a  specific  set  of  loads 
and  operating  conditions. 

The  terms  of  the  ABD  matrix  are  functions  of  the  material 
properties  and  the  ply  angle  0.  If  the  properties  are  kept 
constant  and  within  typical  values  of  most  fiber  resin  types 
of  composite,  it  is  possible  to  plot  the  values  of  Qij  versus 
0  and  get  a  set  of  seven  very  useful  graphs.  These  will  permit 
an  initial  appreciation  of  such  effects,  and  how  they  can  be 
combined  to  meet  a  specific  structural  requirement. 

Figures  5  to  11  show  how  the  elements  of  the  "material 
stiffness"  matrix  ABD  can  be  tailored  by  changing  the  ply 
angle,  taking  full  advantage  of  this  characteristic  behavior 
of  laminated  composite  material . 
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In  figure  11  all  Q's  are  superimposed,  on  the  same  plot, 
to  give  a  better  perspective  of  the  trends  in  doing  this  kind 
of  layer  design  to  fit  specific  structural  constraints.  When 
A16  and  A2g  are  non-zero  there  is  stretching-shearing  coupling, 
while  non-zero  B16  means  twisting-stretching  coupling,  and  non- 
zero B26  indicates  bending-shearing  coupling;  bending-twisting 
coupling  is  originated  in  non-zero  values  of  D16  and  D26  • 

In  normal  design  procedures,  these  terms  are  made  zeroes 
by  selecting  appropriate  stacking  sequence;  cross  ply  (0°/90° 
plies  combination)  or  angle  ply  laminates  (-0/+G  plies 
combination).  Other  important  factor  to  be  considered  is  the 
symmetry  with  respect  to  the  mid-surface  plane. 
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The  stress  strain  relations  or  constitutive  relations  ( ABD 
matrix),  combined  with  the  proper  equations  of  equilibrium  and 
the  strain-displacement  relations  form  the  basis  for  analyzing 
beams,  plates  or  shells. 

These  basic  relations  lead  to  the  motivation  and  reasoning 
behind  the  model  tailoring  to  fit  a  design  requirement. 
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Ill .      ROTOR   MODELING 

A.  TYPICAL   CURRENT   ROTOR 
MODELING 

Current  rotor  design  is  a  compromise  between  the  different 
requirements  in  hover  and  forward  flight,  and  the  analytical 
models  reflect  that  tendency  or  limitation.  For  instance,  the 
optimum  twist  for  hover  is  not  the  same  in  forward  flight  for 
tilt  rotors  as  well  as  for  helicopters.  Indeed,  very  often, 
the  designers  avoid  extension-twist  coupling  with  laminated 
composite  materials  because  of  certain  undesirable  effects. 
The  result  is  that  one  of  the  two  flight  modes  will  have  less 
than  optimal  propulsive  efficiency,  depending  on  design  phase 
decisions  driven  by  requirements  and/or  aircraft  mission 
profiles  . 

With  the  tilt  rotor  concept  of  mixing  fixed  wing  and 
rotorcraft  technology,  these  design  trends  were  brought  to 
surface  once  more. 

B .  PROPOSED   SOLUTION 

Mxon[1987]  proposed  a  solution  for  this  sort  of  dilemma. 
The  model  used  an  extension-twist  coupled  rotor  blade  with 
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100%  RPMD  in  hover  and  80%  RPM  in  forward  flight.  These  are 
reasonable  ratios  between  hover  and  forward  flight  for  a  tilt 
rotor  type  aircraft.  The  difference  in  RPM  provides  a  change 
in  centrifugal  force,  inducing  a  change  in  the  twist  angle 
distribution,  which  can  be  used  to  optimize  the  rotor  for  both 
regimes . 

The  use  of  composite  materials  makes  it  possible  to  design 
a  rotor  blade  that  takes  advantage  of  the  changing  loads, 
rendering  the  required  twist  distribution  optimum  in  each 
flight  mode.  At  this  point  one  question  arises: 

What   is   the  magnitude  of  twist  available  from  a  particular 
geometry   and  material    system? 

The  answer  to  that  question  lies  in  a  method  of 
prediction,  with  a  model  that  behaves  as  close  to  the  actual 
blade  as  possible  yet  inexpensive  enough  to  permit  a  wide 
spectrum  parametric  study,  providing  a  source  to  evaluate 
several  geometry  and  material  systems.  Certainly  such  modeling 
should  be  verified  by  experimental  data. 

This  research  is  aimed  to  address  the  analytical  part  of 
the  above  question  and  provide  design  tools  to  rotor  blade 
systems . 


5  Reference  value  for  100%  RPM  is  217  rad/s,   in  that 
scale  of  rotor  blade. 
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The  torsional  moments  that  affect  a  rotating  blade  can 
modify  substantially  the  twist  distribution. 

In  this  numeric  model,  the  pitching  moments  generated  by 
aerodynamic  loading  are  not  taken  into  account,  the  emphasis 
being  on  the  inertia  and  extension-twist  coupling  as  sources 
of  torsional  moments. 

An  assumed  lift  distribution  acts  on  the  model,  but 
compared  with  the  inertia  forces,  its  effects  are  negligible 
on  the  rotating  blade.  The  inertia  forces  acting  on  the 
rotating  blade  produce  extension  which  causes  twist  by 
structural  coupling.  However,  the  same  inertia  forces  produce 
a  centrifugal  flattening  [Prouty,  1986]  effect  called  tennis 
racket  effect  (because  of  the  tendency  of  a  tennis  racket  to 
align  its  plane  with  the  plane  of  rotation  as  it  is  swung  in 
an  arc ) . 

The  forces  acting  in  a  helicopter  blade  can  be  seen  in 
Figure  12,  and  it  may  be  noted  that  in  this  case  the  effect 
is  reduced  with  balancing  weights  (Chinese  weights). 

C.      GIFTS   CAPABILITIES 

Gifts  capabilities  include  generation  of  any  type  of 
structural  model,  with  the  following  elements: 
.  one-dimension  rods  or  beams 
.  two-dimension  plates 
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three-dimension  shells,   solid  and  complex  stiffened 

shells . 

The  elements  can  be  selected  from  a  library  of  options. 

The  materials  can  be  created  by  the  user  or  selected  from 
a  library  of  standard  definitions.  The  program  allows  the  user 
to  define  isotropic  materials  and  laminated  composite 
materials,  that  can  have  up  to  a  maximum  of  hundred  layers. 

Loads  can  be  concentrated  or  distributed,  static  as  well 
as  function  of  time. 

Static  and  dynamic  analysis  are  performed  on  the  models 
generated.  The  dynamic  analysis  provides  free  vibration 
response  and  forced  response  (transient  computations),  and 
model  superposition. 

In  using  GIFTS  as  the  finite  element  program  to  model  the 
blade,  the  effects  of  the  centrifugal  force  are  taken  into 
account  with  the  forces  acting  on  the  lumped  masses  as  shown 
in  Figure  12,  radially  from  the  center  of  rotation.  In  using 
other  finite  element  programs  this  feature  must  be  verified 
for  correct  usage  and  application.  The  tennis  racket  effect 
can  affect  the  twist  angle  to  a  large  extent.  If  the 
centrifugal  forces  act  in  parallel  lines  instead  of  passing 
through  the  center  of  rotation,  there  will  be  no  centrifugal 
flattening  effect  and  results  in  larger  twist  angles. 

The  basic  rotor  model  selected  to  apply  the  proposed 
method  is  the  one  analyzed  by  Hodges  f  1 987  ]  ,  this  model  was 
analyzed  using  two  different  methods,  namely  by  finite  element 

30 


Chinese 
Weigh! 


Centrifugal 
Force  Acting 
on  Blade  Mass  ■ 
Element 


Center  ol  Rotation 


Blade  Feathering  Axis 


Centrifugal 
Force  Acting  on 
Chinese  Weight 


(MSC/NASTRAN)  and  a 
coupled-beam  model 
[Rehfield,  1985].  The 
rotor  geometry 
presented  by 
Hodges[1987l  was 
modeled  in  GIFTS  using 
triangular  TB3  element 
(  "  flat  triangular 
pi ate  o  r  she 1 1 
element,  including 
membrane,  out  of  plane 
and  fictitious  in- 
plane  bending 
stiffness. "[Batoz, 
19801 ) .  The  results 
indicated    a    wide 

Figure  12  Twisting      Moments      Due      to     disparity  of  results. 

Centrifugal    Forces    [Prouty, 1986 ] .  However,   the  present 

model   appears   to   be 

within   the   assumption  of   linear   and   small   displacements 

theory . 

The  models  in  Figure  13  have  a  notation  (M120  /  Ml-20) 

that   is   described   in   the   next   section.   A   twist   angle 

distribution  discrepancy  could  be  explained  by  the  tennis 

racket  effect  probably  not  present  in  the  model  described  in 


End  View 
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Hodges f 1987 ] ,  but  the  tip  vertical  displacement  due  to  lift 
and  weight  with  no  rotation  could  not  be  explained  by  this 
argument.  It  may  be  observed  that  the  displacements  quoted  in 
the  referenced  work  are  about  25%  of  the  length,  which  appears 
to  violate  the  linear  theory  that  was  used. 


BLADE  DEFLECTION  DUE  TO  LIFT  AND  WEIGHT 


COMPARISON  WITH  NA.5A  TM  83024 


Figure  13  Comparison  between  GIFTS  model  and  NASA  TM  89024 
resul ts. 


This  discrepancy  may  be  due  to  a  possibility  that  the 
models  compared  are  not  exactly  the  same,  although  every 
effort   was   made   to   take   into   account   all   the   pertinent 
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geometric  and  material  data  reported.  However  this  discrepancy 
was  the  driving  factor  that  raised  some  questions,  such  as: 
What  are  the  possibil ities  to  model  different  structural 
configurations  starting  with  the  same  number  of  layers, 
same  lay-up  design,  same  thickness  and  material 
properties?  (perhaps  by  inadvertent  oversight  by  the 
designer,  analyst  using  canned  programs,  or  manufacturing 
engineer. ) 

The  answer  to  this  was  based  on  the  following  assumptions: 
.  The  structure  is  a  D  shape  spar. 

It  is  built  using  three  shells,  top,  bottom  and  rear 
she J  1  surfaces  . 

.  The  top  and  bottom  meet  at  the  leading  edge  and  the  rear- 
shell  closes  the  box. 

.  Once  established,  the  internal  lay-up  in  the  stacking 
sequence  remains  constant . 

With  those  assumptions  a  laminate  to  constitute  one  shell 
surface  may  be  defined,  using  only  the  Laminate  Principal 
Orientation  (  LPO)  angle.  For  example,  in  the  following 
laminate  [  +20°/-  70°/  +  20°/-  70°/-  70°/  +  20°  J  ,  the  top  layer  angle 
is  denoted  +20°  and  is  a  six  layer  asymmetric  laminate. 

Rotating  the  whole  laminate  40°  clockwise  gives 
[-20o/  +  70o/-20o/  +  70f7  +  70o/-20°]  ,  whose  Laminate  Principal 
Orientation  angle  is  denoted  as  -20°. 

This,  laminate  has  intentionally  designed  asymmetry  to  get 
the  extension-twisl  coupling.  Associated  with  this  laminate 
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is  a  normal  vector  whose  orientation  is  seen  in  Figure  14.  The 
D-spar  then  is  built  using  three  such  laminates,  with  each 
laminate  mounted  with  the  normal  vector  facing  inward  or 
outward.  This  method  of  construction  results  in  six  different 
configurations  of  the  rotor.  These  configurations  and  other 
variations  within  each  one  will  be  developed  in  the  next 
section . 

D.       VARIOUS    ROTOR   DESIGNS    OR 
CONFIGURATIONS 

The  six  different  configurations  are  designated  Ml  to  M6 
and  the  LPO6 ,  varying  from  -90°  to  +90°  in  steps  of  20°, 
provides  the  suffix  to  complete  the  model  designation  Ml-907, 
M340,  M6-60  etc.  In  a  total  of  11  LPO  in  each  model. 

Different  configurations  may  be  visualized  with  the  help 
of  the  scheme  shown  in  Figure  14. 

Model  I  and  IV  can  be  built  with  only  one  laminate  wrapped 
around,  one  inverse  of  the  other. 

In  the  program  GIFTS,  the  outward  vector  for  each  surface 
is  defined  during  the  "GRID"  definition,  e.g.  in  a  four  sided 
grid,  each  shell  is  defined  within  a  grid.  The  local  X  axis 


6  Laminate  Principal  Orientation  Angle,  referred  to  the 
top  layer  of  the  laminate  used  to  construct  that  particular 
model . 

7  Ml-90  meaning  model  1  with  LPO  =  -90°;  M340  meaning 
model  3  with  LPO  =  40°. 
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is  oriented  parallel  to  the  first  line  of  the  grid,  the  local 
Z  axis  is  the  cross  product  of  the  first  and  second  "line"  in 
the  grid  definition,  while  the  local  Y  axis  is  the  cross 
product  of  the  Z  and  X  axes,  as  is  usual  in  cartesian  systems. 
The  outward  vector  is  oriented  in  a  direction  defined  by 
the  local  Z  axis.  Layer  angles  are  positive  in  the  counter- 
clockwise direction  in  the  XY  plane  for  shell  elements,  as  is 
the  convention  for  laminated  composite  materials. 
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END  VIEW  OF  D  SPAR  (  3  SHELLS,' 


MODEL' 1   MODEL  2 


MODEL  1 


MODEL  5 
NORMAL  VECTOR  ASSOCIATED  TO 
LAMINATE   STACKING   SEQUENCE 


MODEL  3 


-B> 


MODEL  6 


Figure  14  Normal  to  surface  vector  associated  with  the 
layer  sequence  of  the  basic  laminate  and  its  position  in 
relation    to    end    view   of    the    blade. 


It  may  be  observed  that,  when  an  asymmetry  is  present  in 
the   basic   laminate   (material   definition),   the   structural 
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response  is  highly  dependent  on  how  the  grid  and  consequently 
the  lay-up  is  defined.  In  the  case  of  a  symmetric  laminate  or 
isotropic  material,  this  fact  is  never  brought  to  surface 
because  of  its  irrelevancy.  However,  dealing  with  asymmetric 
laminated  composite  materials,  the  analyst  must  always  have 
this  in  mind.  It  may  be  noted  that,  in  most  cases  the  laminate 
is  chosen  to  be  symmetric. 

Table  1  gives  in  a  compact  form  the  description  of  various 
laminate  constructions  described.  Included  in  that  table  are 
the  line  names  actually  employed  in  the  input  file  M$**.SRC 
[  Annex  1]  ($  =  1,2,3,4,5,6  models;  **  =  top  layer  angle  as 
defined  previously,  in  the  range  -90  to  +90),  to  develop  the 
models  generated  in  GIFTS.  By  modifying  two  areas  in  this 
input  file,  six  models  are  generated  that  are  considered  under 
the  present  study,  this  also  provides  variation  of  the  layer 
angle  within  each  model.  The  actual  modifications  employed  in 
the  input  file  are  given  in  LAYUP.AUX  f  Annex  2].  For  each 
model  (one  of  the  six  combinations  obtainable  by  the  change 
in  the  laminate  associated  vector,  as  given  in  Figure  14), 
the  whole  laminate  was  rotated  in  increments  of  10°, 
referenced  to  the  top  layer  angle  ( LPO )  .  For  each  of  these 
models,  the  RPM ,  and  consequently  the  centrifugal  force,  was 
varied  in  the  range  (130  rad/s  to  300  rad/s). 

The  model  geometric  characteristics  are  the  same  as  in 
Hodges [ 1987  ]  ,  a  graphite  epoxy  composite  rotor  blade 
constructed  out  of  Hercules  IM6  fiber  with  Ciba-Geigy  R6376 
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resin.  The  cured  ply  thickness  is  assumed  to  be  0.0055  in.  The 
D-spar  has  a  35.23  in.  radius,  and  begins  at  station  5.32  in., 
which  has  all  degrees  of  freedom  suppressed  to  emulate  a  rigid 
rotor  connection  to  the  hub.  The  rotor  has  a  constant  cross- 
section,  defined  by  13  nodes  and  is  divided  into  21  span-wise 
(3.0  in  apart)  stations.  Models  with  finer  meshes  were  used 
in  the  early  stages  to  check  for  convergence.  The  mesh 
selected  seemed  to  be  a  reasonable  compromise  between  the 
computational  time  and  accuracy. 

The  orthotropic  material  properties  of  the  lamina  are 
listed  below: 

En,psi 23.1  x  106 

E22,psi 1.4  x  106 

v12 0.3  38 

G12,psi 0.73  x  106 
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TABLE  1  -  MODEL  DEFINITION 


Example 

Vector 

+  20   i 

I 

-70 

+  20 

-70 

-70 

+  20 

POSSIBLE 

MODEI 

jS 

GRIDS,  POINTS  AND  LINES 
L45      4     L14  1 


L51 


L512 


G3 

L411 

Gl 

L18 

G2 

L512 

12    L1112   11     L811         8      L128      12 
(not  to  scale,  e.g.  line  L512  is  30.0  in  length, 

L45  is  .286  in,  and  L14  is  about  1.15  in) 

*  =  vector  inward   ;   o  =  vector  outward. 

RIGHT  HAND  RULE  NORMAL  VECTOR  TO  SURFACE 

GRID  DEFINITION 


Ml$$ 

VECTOR 

* 

* 

* 

GRID 
DEF. 

L512,L1112, 
L411 ,L45// 

L411 ,L811,L18, 
L14// 

L18,L128,L512, 
L51// 

o 

M2$$ 

VECTOR 

* 

* 

GRID 
DEF. 

L512.L1112, 

L411 ,L45// 

L411 ,L811 ,L18, 
L14// 

L512,L128,L18, 
L51// 

M3$$ 

VECTOR 

* 

o 

o 

GRID 
DEF. 

L512,L1112, 
L411 ,L45// 

L18.L811 ,L411 , 
L14// 

L512,L128,L18, 
L51// 

M4$$ 

VECTOR 

o 

o 

o 

GRID 
DEF. 

L411 ,L1112, 
L512.L45// 

L18.L811 ,L411 , 
L14// 

L512,L128,L18, 
L51// 

M5$$ 

VECTOR 

o 

o 

* 

GRID 
DEF. 

L411 ,L1112, 
L512,L45// 

L18,L811 ,L411 , 
L14// 

L18,L128,L512, 
L51// 

M6$$ 

VECTOR 

o 

* 

* 

GRID 
DEF. 

L411 ,L1112, 
L512,L45// 

L411 ,L811 ,L18, 
L14// 

1 1 C T-I-. — T ..J...-X  —  -iJ 

L18,L128,L512, 
L51// 
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IV.      STATIC   AND   DYNAMIC 

ANALYS I S 

The  six  models  described  in  chapter  III  are  basically 
subjected  to  three  main  types  of  independent  static  loads. 
Further,  a  modal  analysis  is  also  performed  to  obtain  the 
first  four  fundamental  modes  of  vibration  frequencies  and  mode 
shapes . 

The  first  load  applied  is  a  torque  at  the  tip,  not 
combined  with  any  other  load. 

The  second  load  case  consists  of  combined  lift  and  blade 
weight  at  zero  RPM ,  i.e.  with  no  centrifugal  forces. 

The  third  load  case  is  the  centrifugal  force  due  to 
rotation  of  the  blade.  The  lift  and  blade  weight  are  present 
and  the  RPM  varies  from  130  rad/s  to  300  rad/s. 

A.      TORSIONAL   LOAD 

A  moment  of  26.4  in-lb  is  applied  distributed  among  the 
nodes  at  the  tip  as  concentrated  moments. 

This  load  case  reveals  the  influence  of  the  laminate 
principal  orientation  angle  ( LPO  )  ,  and  the  torsional  stiffness 
for  each  model.  The  different  response  of  the  six  different 
models  also  is  presented  for  comparison  purposes. 
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The  results  are  presented  in  Figure  15  for  the  six  models 
for    LPO's    of    +20°    and    -20°    in    each    model. 

As  much  as  40%  variation  in  the  response  may  be  observed 
in    the    comparisons. 


ANGULAR     D  I  SPLACEMENT     COMPAR  I  SON 


TORQUE    APPL I  CO  AT  TIP  OF  BLADE 
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MODEL    6 


Figure  15  Angular    Displacement     at     the     Tip    Due     to    Pure 
Torque. 


B.     LIFT   AND   BLADE   WEIGHT   LOADS 

This  loading  is  designed  to  get  the  combined  effect  of  the 
lift  and  blade  weight  given  by  Hodges T 1 987  ]  ,  expressed  in 
terms  of  radial  distance  of  the  center  as  follows: 

(17)     1  =  (0.02222  lb/in)  x  r  -  0.0123  lb/in 
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This  load  is  applied  vertically  in  two  spanwise  rows  at 
the  leading  and  trailing  edges,  in  a  proportion  that  produces 
zero  moments  at  the  quarter  chord.  This  load  distribution  is 
more  homogeneous  and  induces  less  cross-section  distortion, 
as  opposed  to  the  application  along  a  line  on  upper  surface 
as  in  beam  type  elements.  This  distributed  line  loading 
produces  small  torsional  displacement  (10"4  degree). 

The  results  are  presented  in  Figure  16  for  the  six  models, 
with  2  values  of  LPO  (+20°  and  -20°)  for  each  model. 
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Figure  16  Displacement    at    the    Tip    Due    to    Vertical    Loads 
Along    the    Blade. 
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C.      CENTR.I  FUGAL   LOAD 

Next  loading  case  investigated  reveals  the  effect  of 
varying  centrifugal  load,  combined  with  changing  the  LPO  angle 
within  each  model.  The  response  yields  twist  angle  induced  by 
extension-twist  coupling  present  in  the  structure  due  to 
asymmetric  stacking  sequence. 

The  LPO,  referred  to  the  top  layer  angle,  varies  in  the 
range  of  +/-  90°,  the  internal  lay-up  in  the  stacking  sequence 
remaining  constant,  as  previously  defined. 

After  some  experiments  to  localize  inflection  points  in 
the  twist  angle  variation,  the  LPO  to  be  investigated  were 
chosen  in  the  range  of  +/-  90°,  giving  a  total  of  11  laminate 
principal  orientation  angle  for  each  of  the  six  models. 


Y7>97\ 


■••  • 


vaa_ 


///A 


~, 


//AS 


V>V^ 


V 


yy 


/2a 


Y/s 


ga 


///A Y/A 


Y/A/A Y^U 


V/s 


w, 


//A/a 


Y/A/A 


'/■ 


/^ 


^L 


MOOEl.  2       M0D6L  3 


Figure  17  Overall    Range    of    the    Twist   Angle    within    the   Six 
Model s . 
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The  loads  due  to  lift  and  blade  weight  are  present  but 
their  effects  are  observed  to  be  negligible  compared  with 
those  due  to  centrifugal  loads. 

The  RPM  is  varied  from  130  rad/s  to  300  rad/s  in 
increments  of  10  rad/s,  in  a  total  of  19  load  cases.  The 
combination  of  load  cases  and  LPO  give  a  total  of  201  points 
of  twist  angle  for  each  of  the  six  models.  These  points  are 
presented  in  six  3-D  graphs  to  better  appreciate  the  choices 
offered  in  this  design  approach.  Those  graphs  present  the 
trends  of  the  twist  angle  with  RPM  and  LPO  variations.  The 
twist  angle  range  for  each  model  is  in  Figure  17,  while  the 
overall  parametric  information  are  presented  in  Figures  18 
t  h  rough  23  .. 


Figure  18  Model  One,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),    Referred    to    Top   Layer  Angle. 
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Figure  19  Model  Two,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
( LPO ) ,    Referred    to    Top   Layer  Angle. 
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Figure  20  Model  Three,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
( LPO ) ,  Referred    to    Top    Layer   Angle. 
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Figure  21  Model  Four,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
( LPO ) ,     Referred    to    Top    Layer   Angle. 
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Figure  22  Model  Five,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),    Referred    to    Top   Layer  Angle. 
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Figure  23  Model  Six,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),    Referred    to    Top   Layer   Angle. 
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3D  . 


MODAL   ANALYSIS 


To  observe  the  effects  of  each  configuration  (Model)  and 
laminate  principal  orientation  on  the  first  four  fundamental 
frequencies,  a  free  vibration  analysis  is  performed. 

The  results  of  the  modal8  analysis  are  depicted  in  Figures 
24  and  25  for  the  six  models  and  for  LPO's  of  -20°  and  +20° 
respectively . 
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\s'  .'  I   1  Fvi 


MODE^  2       MODEL  3        MODEL  <* 
FLATWISE  AND  CHORDWlSE  MODES 


IV  VI   -I  cw 


W?l       2  F* 


ES  2  c* 


Figure  24  Frequency    Variation    among   Different    Models 


1  FW  means  First  Flapwise;  1  CW  means  First  Chordwise 
and  so  on. 
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It  may  be  noted  that,  even  for  simple  models,  this 
solution  procedure  is  a  time  consuming  task,  only  recommended 
for  final  or  quasi  final  configurations  in  the  design  process, 
or  using  a  very  crude  mesh  lay-out. 

It  can  be  seen  that,  for  each  vibration  mode,  the 
frequency  remains  practically  constant.  This  reveals  a  low 
sensitivity  of  the  frequency  to  changes  in  model 
configurations,  considering  the  first  four  fundamental  modes 
and  same  LPO ' s . 
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Figure  25  Frequency    Variation    among   Different    Models. 
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V  .      ANALYSIS   OF   RESULTS 


In  this  chapter  the  results  obtained  for  each  model  will 
be  presented  and  several  observations  are  inferred. 

The  present  research  was  conducted  to  get  an  insight  of 
the  trends  implicated  into  the  two  questions  posed  in  Chapter 
III,  sections  B  and  C.  This  study  also  provides  some  insight 
into  answering  the  questions  posed  in  sections  B  and  C  of 
Chapter  III. 


What   is   the  magnitude  of  twist  available  from  a  particular 
geometry   and  material    system? 


What  are  the  possibilities  to  model  different  structural 
configurations  starting  with  the  same  number  of  layers, 
same  lay-up  design,  same  thickness  and  material 
properties ? 


The  first  question  is  a  design  approach  to  structural 
analysis.  That  approach  will  be  elaborated  in  this  chapter. 
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The  second  question  opened  a  vast  field  for  research 
activity  that  will  be  addressed  in  next  Chapter. 

It  can  be  seen  in  Figure  14  that  the  models  1  through  6 
are  obtained  simply  by  changing  one  of  the  normal  vectors 
associated  with  the  laminate,  (in  one  of  the  three  faces  of 
the  spar).  The  effects  of  these  changes  in  the  structural 
behavior  of  the  blade  may  be  divided  into  two  main  categories, 
static  and  dynamic  behavior. 

A.      DYNAMIC   BEHAVIOR, 

The  free  vibration  characteristics  are  shown  in  Figures 
24  and  25.  Figure  24  shows  the  frequencies  for  LPO  of  -20°  and 
Figure  25  presents  for  LPO  of  20°.  The  first  four  fundamental 
modes  of  vibration  have  very  little  sensitivity  to  the  model 
type,  and  may  be  observed  that  the  frequencies  do  not 
change,   for  the  same  or  symmetric  LPO  (positive/negative). 

This  behavior  can  be  considered  benign,  or  in  other  words 
auspicious,  in  a  connotation  that,  once  the  designer  arrived 
at  the  desired  frequency  values  for  the  RPM  range  considered, 
the  LPO  or  the  model  can  be  changed  until  a  desired  twist- 
extension  coupling  is  achieved. 

These  iterations  can  be  done  without  looking  very  close 
to  resonance  problems,  which  must  be  checked  anyway  when  the 
final  model  is  to  be  selected. 


This  assumption  is  made  based  on  the  validity  of 
Yntema[1955]  approach  for  these  models.  It  is  possible  that 
the  change  in  the  stiffness  matrix  due  to  RPM  affects  the 
rotor  built  of  laminated  composite  materials  in  a  different 
way  than  isotropic  rotors.  Further  research  needs  to  be  done 
along  this  direction  to  obtain  comparable  results  for 
composite  rotors. 

In  Figures  24  or  25,  each  column  corresponds  to  one  mode 
of  vibration  at  zero  RPM  and  is  a  point  in  the  3-D  plot  for 
that  specific  mode  (flatwise,  chordwise  or  pitchwise9).  The 
RPM  and  LPO  may  be  varied  as  in  the  static  case  yielding 
carpet  plots.  The  resulting  plots  would  validate  or  not  the 
assumption  of  Yntema's  values  for  composite  blades. 

B.      STATIC   BEHAVIOR, 

The  static  behavior  is  presented  in  Figures  18  through  23. 
The  parameters  considered  are  RPM  and  LPO  as  inputs  and  the 
twist-extension  coupling  manifested  through  the  twist  angle 
at  the  tip,  as  the  output.  This  angle  will  be  referred  as  TTA 
(tip  twist  angle)  in  the  subsequent  discussion. 

The  TTA  can  be  used  as  a  measure  of  the  twist  angle  ratio 
along  the  blade  span.  The  twist  angle  is  not  a  linear  function 
of  the  radius  but  may  be  represented  by  a  second  degree 


9  Pitchwise  mode  is  a  higher  mode  in  this  case  and  is  not 
considered  in  this  report. 
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polynomial,  such  that  the  TTA  is  a  good  parametric  factor  to 
infer  the  twist  angle  magnitude  along  the  blade. 

As  mentioned  before,  the  six  models  are  obtained  by 
changing  the  orientation  of  the  normal  vector  associated  to 
the  laminate,  in  the  three  shells  that  constitute  the  spar, 
one  at  a  time  as  shown  in  Figure  14. 

In  that  figure  it  can  be  verified  that  Model  1  and  Model 
4  can  be  built  with  only  one  unbroken  piece  of  laminate, 
wrapped  around  a  mold.  Further,  one  model  is  just  the  inverse 
of  the  other,  in  other  words  Model  1  is  the  inside  out  of 
Mode]  4. 

The  models  2  and  5  can  also  be  considered  as  some  sort  of 
inverse,  in  a  broader  sense;  similar  reasoning  may  be  extended 
to  models  3  and  6. 

With  the  above  considerations,  the  six  models  can  be 
analyzed  in  pairs,  which  is  very  convenient  from  the  point  of 
view  of  getting  the  information  in  a  simple  and  concise  form, 
so  that  the  designer'  can  quickly  evaluate  the  model  ,  keeping 
in  view  the  design  requirements. 

The  second  important  consequence  of  this  inverse  concept 
is  that  ,  Mode  J  1  and  Mode]  4  can  supply  a  spectrum  of  twist 
angle  requirements,  for  a  given  loading  condition. 

The  importance  of  this  reside  in  the  manufacturing 
process,  in  that  both  models  may  be  built  by  filament  wound 
technique.  This  process  is  less  expensive  arid  more  reliable 


as  it  eliminates  the  filament  discontinuity  in  the  shell's 
junction,  and  requires  less  accurate  quality  control. 

The  convention  used  in  the  modeling  process  has  as  a 
aftermath  that  positive  angles  have  active  aerodynamic 
stabilizing  effect,  acting  as  a  propeller  governor  device. 

1 .  Models  1  and  4 

Figures  18  and  24  exhibit  the  behavior  of  these  two 
model s . 

The  models,  which  were  referred  as  inverse  to  each  other, 
also  show  the  response  characteristics  and  behave  as  expected. 

This  may  be  noted  by  observing  the  displacement  for  a 
given  RPM  and  Layer  Principal  Orientation  Angle  in  Figure  26. 

2.  Models  2  and  5 

Those  models  have  the  same  inverse  relationship  in 
the  TTA  value.  The  main  difference  is  that  these  two  models 
are  stiffer  than  the  other  four,  as  shown  in  Figure  17.  Figure 
2  7  shows  these  two  models  drawn  together. 

3.  Models  3  and  6 

The  inverse  relationship  also  holds  here  and  the 
stiffness  values  are  higher  than  that  found  in  models  1  and 
4  and  lower  than  models  2  and  5  (Figure  17).  See  Figure  28 
for  this  compai  is on. 
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MODEL  I 


&* 


Figure  26  Models  1  and  4  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  ( LPO )  angle, 
output    is    Tip    Twist    Angle    ( TTA ) . 
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Figure  27  Models  2  and  5  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  ( LPO )  angle, 
output    is    Tip    Twist    Angle    ( TTA ) . 
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Figure  28  Models  3  and  6  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  ( LPO )  angle, 
output    is    Tip    Twist    Angle    (TTA). 
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VI  .      CONCLUSIONS   AND   SCOI>E 
FOR   FUTURE   RESEARCH 


A  . 


GONCLUS IONS 


This  report  presents  a  detailed  analysis  of  a  typical 
rotor  blade,  the  focus  being  the  static  behavior  in  the 
presence  of  extension-twist  coupling  due  to  asymmetric 
stacking  sequence  of  laminas. 

Different  extension-twist  coupling  effects  are  obtained, 
by  changing  the  ply  orientation  with  respect  to  the  body  axes, 
while  keeping  the  stacking  sequence,  thickness  and  internal 
lay-up  sequence  as  invariants. 

These  effects  are  presented  in  terms  of  carpet  plots.  The 
extension-twist  coupling  is  measured  through  the  angular 
displacement  at  the  tip  of  the  blade  as  function  of  RPM  and 
Laminate  Principal  Orientation  (  LPO  )  angle  variations. 

Free  vibration  analysis  is  performed  for  two  LPO's  of  the 
six  different  models.  The  first  four  fundamental  frequencies 
show  small  sensitivity  to  these  specific  structural  changes. 

The  static  analysis  revealed  a  sort  of  "antisymmetric" 
behavior  within  the  models,  allowing  them  to  be  grouped  in 
pairs.  Such  behavior  can  be  of  use  in  the  design  process,  to 
trim  structural  response  to  prescribed  loads. 
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A  structural  analyst  should  be  very  careful  in  dealing 
with  composite  laminated  structures,  especially  if  the  analyst 
is  experienced  only  in  structural  analysis  of  isotropic 
materials.  It  was  noted  that  it  is  very  easy  to  be  mistaken, 
mainly  if  there  are  asymmetries  present  within  the  laminate. 
It  is  very  important  to  know  how  the  local  coordinates  systems 
are  established  in  the  finite  element  program  being  used. 
These  axes  will  determine  the  relative  position  of  the 
laminate  with  respect  to  the  structure.  If  asymmetries  are 
present,  it  is  possible  to  get  different  structural  response 
and  have  different  structural  configurations,  all  starting 
with  the  same  basic  laminate. 

B  .      SCOPE   EOR   FUTURE   RESEARCH 

Future  research  may  be  beneficial  into  two  main  areas. 

The  first  one  being  to  validate  these  models  with 
experimental  data,  providing  the  necessary  confidence  in 
further  utilization  of  the  method  envisaged. 

Once  the  validation  is  done,  the  method  can  be  corrected 
and  trimmed  for  operational  use. 

Then  a  field  that  calls  for  attention  is  a  dynamic 
characteristics  mapping,  that  is,  with  the  RPM  and  LPO  as 
parameters  obtain  the  fundamental  frequencies  as  output. 

This  would  be  a  3-D  composite  material  extension  of 
Yntema's  [1955]  work  for  beams  of  isotropic  material. 
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ANNEX   1      —      M120.SRC 

$         DEFINES  BLADE  MODEL  WITH  NO  ELIPSE  PROJECTION 
$         (  LOST  3.976  DIGITS  IN  DECOM  ) 

*  (  ~  20  MIN  TO  RUN  2  LDCA5ES  ) 

$         NO  ROOT  DEFINED  JUST  SUPRESSED  DOF  OF  FIRST  RIB 

$ 

KPOINT 

1/ .65,5.23,0 

2/. 455, 5. 23, . 1 103 

3/0,5.23, . 1545 

4/- .49,5.23, . 1423 

5/- .49,5.23, - . 1433 

6/0,5.23, - . 1545 

7/ .455,5.23, - . 1 103 

8/ .65,35.23,0 

9/ .455,35.22, . 1103 

10/0,35.23, . 1545 

1  1/- .49,35.23,  . 1433 

12/- .49,35.23, - . 1433 

13/0,35.23, - . 1545 

14/  .455,35.23, - .  1 103// 

$         LINES  GENERATION 

$ 

SLINE 

L34/3,4,2 

L45/4,5,4 

L56/5.6,2 

L18/1 ,8,21 

L41 1/4, 11,21 

L512/5, 12,21 

L101 1/10,11,2 

L  1  1  1  2  /  1  1  ,  1  2  ,  4 

L 1213/ 12, 13,2// 

$ 

*  ARC  GENERATION 
$ 

CARC 

CI  3/1  ,2,2,5 

C6 1 / 6 , 7 , 1 , 5 

CS 10/6,9, 10,5 

C138/13, 14,8,5// 

f 

J         COMPOSITE  LINES  GENERATION 

$ 

COMPL INE 

L14/C13,L3- 

L51/L56,C61 

LSI 1/CS10,L 101 1 

L 126/ L  1213,2136 

7L 


$         BLADE  MATERIAL  AND  ElEMENTS  DEFINITION  (IM6/R6376) 

*  BLADE  MASS  DISTRIBUTION  =  .0123  LB/IN  =  .00038233  SLUG/IN 
$         (COMPOSITE  DENSITVr .001796  SLUG/INA3) 

S         (TUNGSTEN  DENSITY  =.021116  SLUG/INA3) 

$         (MEAN  DENSITY      =  . 004418  SLUG/INA3) 

$         (COMPOSITE  AREA    =.086525  INA2) 

$         (COMPOSITE  VOLUME  =2.5957   INA3) 

$         (TUNGSTEN  AREA     =.01358   INA2   =>  DIA=.1315  IN) 

$         (MODEL  FOCUSED  LENGHT=30  IN;  TOTAL  LENGHT=35.23  IN) 

0RMAT,53 

1/5.07E5,23. 1E6, .338, .73E6, .004418 

8. 12E3, 1 .4E6, 1 .41E4// 

$ 

*  LAYER  DEFINITION  [+20/ -70/+20/ -70/ -70/+20 ] 
$         THICKNESS  =  .0055  IN/LAYER  *  6  =  .033  IN 
COMPTH 

1,6/1,1,1,1,1,1 
.0055, .01375,20 
.0055, .00825,-70 
.0055, .00275,20 
.0055, - .00275,-70 
.0055, - .00825, -70 
.0055, -.01 375, 20// 

*  ELEMENTS 

$  GRIDS  MUST  BE  GENERATED  IN  SAME  ORIENTATION 

$  STABLISH  A  VIEW  POINT  NORMAL  TO  SURFACE  OUTSIDE  (OR  INSIDE) 

X  WITHOUT  CHANGING  THIS  REFERENCE  GENERATE  THE  GRIDS  FOR 

$         ORTHOTROPIC  MATERIALS  (FILAMENT  WOUNDED  THINKING) 

$ 

$         THE  FIRST  SIDE  OF  GRID  DETERMINES  THE  REFERENCE  FOR  MATANG 

$ 

GETY/TB3/1 , 1// 

GRID4 

G1/L41  1  ,L61 1  ,L18,L14 

G2/L1B,L128,L512, L51 

G3/L512,L1 1 12.L41 1 ,L45// 

END 


ANNEX       2  —  LAYUP . AUX 

FILE    TO    BE    USED     IN    BUILDING    THE    SEVERAL     LAYUPS    IN 

EACH    M$**.SRC    FILES 

$  =  LAYUP  NUMBER 

**  =  ANGLE  OF  TOP  LAYER  WITH  RESPECT  TO  NORMAL  VETOR  DIRECTION 

1-  LAYUP  GRID  DEFINITION  --  6  MODELS  -- 

Ml**. SRC 
GRID4 

G1/L411,L811,L1S,L14 
G2/L1G, LI  28, LSI  2, LSI 
G3/L512, LI  1 12 , L41 1 ,L45// 

M2** .SRC 

GRID4 

G1/L41 1 , LB1 1 , L IS, L14 

G2/L512, L12S, LIS, L51 

G3/L512, Li  1 12, L4  1  1  ,L45// 

M3**.SRC 

GRID4 

G1/L1B,LS1 1 ,L41 1  ,L14 

G2/L512, L128,L1S,L51 

G3/L512, LI  1 12, L 411 , L45// 

M4**.SRC 
GRID4 

B1/L1S,L811,L411,L14 
G2/L512, L12S, L1B,L51 
G3/L411 ,L1112,L512,L45// 

M5** .SRC 
GRID4 

Gl/L18,LSil,L4H,L14 
G2/L18, L128,L512,L51 
G3/L41 1 ,L11 12, LSI  2, L45// 

M6** .SRC 
GRID4 

Gl/L411,Lbil,LlB,L14 
G2/L IS, L12S, LSI  2 , L51 
G3/ L  4  i  1  , L 1 1 1 2 , L  5 1 2 , L 4S / / 
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2-  TOP  LAYER  ANGLE  --   IB   ANGLES  -- 

MSOO.SRC 
.0055, .01375,0 
.0055, .00825, -90 
.0055, .00275,0 
.0055, - .00275, -90 
.0055, -.00825, -90 
.0055, - .01375,0// 

*******************♦******#****+*♦********** 

MS10.SRC 

.0055, .01375, 10 

.0055, .00825, -80 

.0055, .00275, 10 

.0055,- .00275,-80 

.0055, - .00825, -80 

.0055,- .01375, 10// 


M$20.SRC 

0055, .01375,20 
. 00825 , -70 
. 00 175 ,20 
- . 00275 , -70 
- .00825, -70 
- .01375,20,  . 


0055 
0055 
0055 
0055 

0055 


***#  +  ♦  +  #  +  *********  +  ♦#  +  *  +  ♦♦*♦-**♦*♦♦♦*♦**♦.♦♦  +  * 

M$30.SRC 
. 0055 , . 01375 , 30 
. 0055, .00825, -60 
.0055, .00275,30 
.0055, - .00275, -60 
.0055, - .00825, -60 
.0055, - .01375,30// 


MS40.SRC 
.0055, .01375,40 
.0055, .00825, -50 
.0055, .00275,40 
.0055, - .00275, -50 
.0055, - .00825, -50 
.0055. - .01375,40// 


M$50.SRC 

.0055, .01375,50 
.0055, .00825, -40 
.0055, .00275,50 
.0055, -.00275, -40 
.0055,- .00825,-40 
.0055, - .01375,50// 

**#**********♦*****♦*********+*******♦****** 

MS60.SRC 
.0055, .01375,60 
.0055, .00825,-30 
.0055, .00275,60 
.0055, -.00275,-30 
.0055. -.00825, -30 
.0055, - .01375,60// 

********##*******♦***+*****♦******#********# 

M$70. SRC 

.0055, .01375,70 

.0055, .00825, -20 

.0055, .00275,70 

.0055, - .00275, -20 

.0055, - .00825, -20 

.0055. - .01375,70// 

**  +  ♦*  +  *****♦***  +  ****  +  *  +  ****  +  +  ****♦*-*■*#♦♦*  +  ♦* 

M*8C  . Er  C 
. 0055 , . 01375 , 80 
.0055, .00625,-10 
.0055, .00275,80 
.0055, - .00275, -10 
.0055, -.00825, -10 
.0055, - .01375,80// 

**  +  ****■***■*•*****+-*♦♦**♦■************♦*♦*♦*■»♦* 

MS90.SRC 

.0055, .01375,90 

.0055, .00825,0 

.0055, .00275,90 

.0055, - .00275,0 

.0055 , - . 00525 , 0 

.0055. - .01375.90 


M$-90.SRC 

.0055, .01375, -90 
.  X>55,  .0C  :-  I"  .  . 
'55  ,  .  00275  ,  -5 


*■*#*♦***-*■+■******  +  **  + 


.  -  . 


55  ,  -  .  C     32     . 
.        55, -.01375, 


M*-BO.SRC 

.0055, .01375,-80 
.0055, .00B25, 10 
.0055, .00275,-80 
.0055, -.002 75, 10 
.0055, -.00825, 10 
.0055,- .01375,-80// 

********************************************* 

MS-70.SRC 

.0055, .01375, -70 

.0055, .00825,20 

.0055, .00275,-70 

.0055,- .00275,20 

.0055, - .00825,20 

.0055, -.01 375, -70// 

♦  **#♦•♦******♦♦♦***#♦***  +  ***♦♦**♦  +  **■♦■)<-**♦■»♦■»•* 

MS-60.SRC 

.0055, .01375,-60 

.0055, .00825,30 

.0055, .00275,-60 

.0055, - .00275,30 

.0055.- .00825,30 

. 0055 , - . 0 1 275 , -60// 

MS -50. SRC 
.0055, .01375,-50 
.0055, .00325,40 
.0055, .00275,-50 
.0055, - .00275,40 
.0055, - .00825, AC 
.0055, - .01575,-50// 

**  +  *♦•  +  ******♦*************  +  ****♦*♦  +  ♦♦♦*♦♦♦.*♦ 

MS -40 . SRC 

.0055, .01375, -40 

.0055, .00825,50 

.0055, .00275, -40 

.0055, - .00275,50 

.0055, - .00825,50 

.0055, - .01375, -40// 


M*-30.SRC 

.0055, .01375 , -50 
.0055, .00525,60 
.  0055  .  00275  .  -  30 
.0055. - .00275,60 
.0055, - .00825 , faC 
.  0055 . - , 01 375. -30/ 


•  ♦♦♦HI-**********************************-****-* 

MS-20.SRC 

.0055, .01375, -20 

.0055, .00825,70 

.0055, .00275,-20 

.0055, -.00275, 70 

.0055, - .00825,70 

.0055, -.01375, -20// 

*******+♦***♦*♦**♦*♦***#*****♦*********♦♦*#* 
MS-10.SRC 

.0055, .01375, -10 
.0055, .00825,80 
.0055, .00275, -10 
.0055, - .00275,60 
.0055, -.00825,80 
.0055, - .01375, -  10// 


7£ 


INITIAL       DISTRIBUTION       LIST 


N°  Copies 


Defense   Technical   Information  Center  2 

Cameron  Station 
Alexandria,  VA  22304-6145 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

Centro  Tecnico  Aeroespacial  8 

IPD  -  Divisao  de  Aeronautica 

Attn:  Chefe  da  PAR 

Sao  Jose  dos  Campos,  SP  BRAZIL 

CEP  12225 


Department   Chairman,   Code   67 
Department  of  Aeronautics  and  Astronautics 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Department     of    Aeronautics 

Attn:  Professor  Ramesh  Kolar,  Code  67Kj 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

Department  of  Mechanical  Engineering 
Attn:  Professor  Gilles  Cantin,  Code  69Ci 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Major  Geraldo  A.  Macedo  Moura 
Centro  Tecnico  Aerospacial 
IPD  -  Divisao  de  Aeronautica 
Sao  Jose  dos  Campos,  SP  BRAZIL 
CEP  12225 

David  Taylor  Research  Center 
Attn:  Dr.  Rem  Jones,  Code  1732 
Bethesda,  MD  20084-5000 

Aerof lightdynamics  Directorate 

Rotorcraft  Dynamics 

Attn:  Dr.  Robert  E.  Ormistron 

MS  215-1 ,  NASA  AMES 

Moffett  Field,  CA  94035 


79 


10.  NAVAIR,  Code  PMA  261 

H-53/H-46  Helicopters  Program  Office 
Attn:  Program  Manager 
1411  Jefferson  Davis  Hwy,  JPI , 
Washington,  D.C.  20361-9310 

11.  Department   of  Mathematics 

Attn:  Professor  D.  Danielson,  Code  53Dd 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 


80 


Thesis 
Ml 839 2 

c.l 


Macado  .lour a 

An  approach   f  or,,dlfsign 
and  analysis   o^composi'ca 
rotor  bli 


